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Abstract 

A two-grid scheme based on mixed finite-element approximations to 
the incompressible Navier-Stokes equations is introduced and analyzed. 
In the first level the standard mixed finite-element approximation over 
a coarse mesh is computed. In the second level the approximation is 
postprocessed by solving a discrete Oseen-type problem on a finer mesh. 
The two-level method is optimal in the sense that, when a suitable value 
of the coarse mesh diameter is chosen, it has the rate of convergence of the 
standard mixed finite-element method over the fine mesh. Alternatively, 
it can be seen as a postprocessed method in which the rate of convergence 
is increased by one unit with respect to the coarse mesh. The analysis 
takes into account the loss of regularity at initial time of the solution of the 
Navier-Stokes equations in absence of nonlocal compatibility conditions. 
Some numerical experiments are shown. 



1 Introduction 

We consider the incompressible Navier-Stokes equations 

u t - vAu + (u • V)li + Vp = /, (1) 
div(u) = 0, 

in a bounded domain ft C M. d (d = 2,3) with a smooth boundary subject to 
homogeneous Dirichlet boundary conditions u — on dCt. In ([I]), u is the 
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velocity field, p the pressure, v > the diffusion coefficient and / a given force 
field. 

In this paper we study the following two-grid mixed finite-element method 
for the spatial discretization of the above equations. First, for the solution {u,p) 
of the fully nonlinear Navier-Stokes equations ([JJ) corresponding to a given initial 
condition 

u(.,0) = uo, (2) 

the mixed finite-clement approximation (uh,Ph) over a coarse mesh of diameter 
H is computed. Then, for any time t > 0, the postprocessed approximation 
{uh,Ph) is obtained as the mixed finite-clement approximation over a finer mesh 
(h < H) to the following steady Oseen-type problem: 

-uAu + (u H (t) • V)fi + Vp = f - f t u H (t)) 

div(u) = Oj 

u = 0, 

In this paper we prove that, in terms of the fine mesh diameter h, this two- 
grid technique is of optimal order in the sense that, for appropriate choices 
of the coarse mesh diameter H, the method has the same rate of convergence 
of standard mixed finite element approximations in the fine mesh. On the 
other hand, for a suitable value of the discretization parameter h, the rate of 
convergence of the postprocessed approximation in terms of H increases by 
one unit the rate of convergence of the coarse standard approximation. The 
improvement in precision is achieved in both the H 1 norm for the velocity 
and the L 2 norm for the pressure in the case of linear, quadratic and cubic 
elements. For other than linear elements the rate of convergence in the L 2 norm 
of the velocity is also increased by one unit. We remark that time evolution is 
performed only at the coarse mesh whereas at the fine grid the time appears 
only as a parameter (see equation (j3j), thus the name of static two-grid method. 

Two-grid or two-level methods are a well established technique for nonlinear 
steady problems, see [33]. In [25], [35] several two- level methods are considered 
to approximate the steady Navier-Stokes equations. They require solving a 
nonlinear system over a coarse mesh and, depending on the algorithm chosen, 
one Stokes problem, one linear Oseen problem or one Newton step over the fine 
mesh. The corresponding algorithms obtain the optimal rate of convergence in 
the fine mesh for appropriate choices of the coarse mesh diameter H. 

In the case of nonlinear evolutionary equations, two-grid techniques have 
been proposed and studied in p], [22], [24], [14]. In these methods, as opposed 
to the method studied in the present paper, time evolution is also performed over 
the fine mesh. The advantage of the method studied in the present paper is that 
since the time integration is only carried out on the coarse mesh, computations 
on the fine grid can be done at selected target time levels where an improved 
approximation is desired, with the corresponding reduction of computing time, 
specially if these target time levels are sufficiently spaced in time. For this 
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reason, although some of the two-grid methods that incorporate the evolution 
in time of the fine mesh approximation are more accurate, the method we present 
can still be more efficient in terms of computational effort for a given error level. 

Two-grid techniques that integrate in time only on the coarse level have 
previously been developed in [TD], [UJ (see also [37]) for spectral methods, and 
later extended to mixed finite-element formulations in [3J, [3], [TDJ. In all these 
works the two grid method is referred to as postprocessed Galerkin method, 
and, instead of ([3J, the approximation (v,h,Ph) is found as an approximation to 
the following Stokes problem 

-vAu + Vp = / - f t u H (t) - (u H (t) ■ V)u H {t) } 

div(u) = Oj 

u = 0, 

This two-grid method will be termed standard postprocessed method, to dif- 
ferentiate it to that studied in the present paper, which will be termed new 
postprocessed method. Both, the standard and the new postprocessed meth- 
ods, have the same rate of convergence. However, as already noted in [12] for 
nonlinear convection-diffusion problems, the new postprocessing technique pro- 
duces more accurate approximations than the standard postprocessed method, 
for moderate to small values of the diffusion parameter v. This will also be the 
case in the numerical experiments in the present paper for moderate values of 
the Reynolds number. 

In the present paper we take into account the loss of regularity suffered by 
the solutions of the Navier-Stokes equations at the initial time in the absence 
of nonlocal compatibility conditions. Thus, for the analysis, we do not assume 
the solution u to have more than second-order spatial derivatives bounded in 
L 2 up to initial time t = 0, since demanding further regularity requires the data 
to satisfy nonlocal compatibility conditions unlikely to be fulfilled in practical 
situations [15], [2D]. Due to the loss of regularity at t = 0, the best error bound 
that we can obtain is 0(H 5 \ log(_ff)|). For this reason we do not analyze higher 
than cubic finite elements. The same limit in the rate of convergence was found 
in 20 for standard mixed finite-element approximations and in |10| . [14] for 
two-grid schemes. 

In practice, any method to numerically solve evolutionary equations needs 
of some time discretization procedure. For brevity reasons, we have preferred to 
present the method in a scmidiscrctc manner without reference to any particular 
time discretization. However, we emphasize that being static, the method we 
present can be applied exactly in the same form, with any time discretization. 
The analysis of fully discrete procedures can be developed along the same lines 
that appear in [IT], [T5] , 

The rest of the paper is as follows. In Section 2 we introduce some prelim- 
inaries and notation. In Section 3 we carry out the error analysis of the new 
method. Finally, some numerical experiments are shown in the last section. 
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2 Preliminaries and notations 

We will assume that fl is a bounded domain in M. d , d = 2,3, of class C ,n , for 
m > 2. When dealing with linear elements (r = 2 below) f2 may also be a 
convex polygonal or polyhedral domain. We consider the Hilbert spaces 

H = {u€ L 2 (0) d | div(w) = 0, u ■ n lon = 0} , 
V = {u e Hl{0) d I div(u) = 0} , 

endowed with the inner product of L 2 (0) d and /Zq(0) , respectively. For I > 
integer and 1 < q < oo, we consider the standard spaces, W l ' 9 (Sl) d , of functions 
with derivatives up to order I in L q (Q), and H l (fl) d = W^' 2 (Q) d . We will 
denote by || • ||; the norm in H l (£l) d , and || ■ ||_; will represent the norm of its 
dual space. We consider also the quotient spaces if'(f2)/R with norm \\p\\h> /r = 
inf{||p + c||, \ceR}. 

We recall the following Sobolev's imbeddings [2J: For q G [1, oo), there exists 
a constant C = C(f2, q) such that 

Hi*' < C\\v\\ w *, 9 , I>I-£>0, <z<co, ueW^(n) d . (5) 
q q a 

For q' = oo, © holds with ± < |. 

The following inf-sup condition is satisfied (see [H]), there exists a constant 
f3 > such that 

inf sup „ ( ,f, V / ) >/?. (6) 

Let II : L 2 (0) d — > H be the L 2 (0) rf projection onto We denote by A 
the Stokes operator on 0: 

A:V(A)cH — > H, A = -UA, V{A) = H 2 (0) d n V. 

We shall assume that u is a strong solution up to time t = T, so that 

||u(t)||i<Afi, Kt)||2<M 2) 0<t<T, (7) 

for some constants Mi and M 2 . We shall also assume that there exists a constant 
Mi such that 

||/||i + ]]/ t ||i + ||/ it ||i <M 2 , 0<t<T. (8) 
Finally, we shall assume that for some k > 2 

L(fc-2)/2j 

sup ||a} fc/2j /|| fc _ 1 _ 2 | fc/2 , + V] sup ||9 t J /|| < +oo, 

0<t<T fc 1 2Lfc/2J 0<t<T "* 23 2 

so that, according to Theorems 2.4 and 2.5 in [T!5], there exist positive constants 
Mfe and such that the following bounds hold: 

||«(t)||fc + ||ut(t)||fc-2 + ||p(t)||ff»-i/R < M k T(tf- k ' 2 , (9) 
f a fc _ 3 ( S )(||u( S )|| 2 + \\u s ( S )\\ 2 _ 2 + ||P(S)||VVE + IIP-(*)II^-s/k) ds < K l ( 10 ) 
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where r(f) = min(i, 1) and a n — e _Q ^" s ^r"(s) for some a > 0. Observe that 
for t < T < oo, we can take r(t) = t and <x„(s) = s". For simplicity, we will 
take these values of r and a n . 

Let 7/i = (Tj- ,<fri)i£i h , h > be a family of partitions of suitable domains 
where /i is the maximum diameter of the elements t' 1 G 7~h, and <^ are the 
mappings of the reference simplex To onto r/ 1 . 

Let r > 2, we consider the finite-element spaces 

S h , r = [xh e C (0 /t ) | xh| T * o ^ e P r - x (To)} c H\Q H ), s° Kr = S h , r nH^0 h ), 

where P r_1 (r ) denotes the space of polynomials of degree at most r — 1 on t . 

We shall denote by (Xh, r , Qh,r-i) the so-called Hood-Taylor element [71 [21] . 
when r > 3, where 

Xh,r = (Sl r ) d , Qft >r -i = S fc , r _i n L 2 (0 h )/R, r > 3, 

and the so-called mini-element [5] when r = 2, where = >S/i,2 H L 2 (0/ l )/R, 
and Xh,2 = {S\ 2 ) d ®V>h- Here, is spanned by the bubble functions b T , r G 7/j, 
defined by & T (x) = (ci + l) d+1 Ai(a;) • • • \d+i{x), if a; € r and elsewhere, where 
Ai(a;), . . . , Xd+i(x) denote the barycentric coordinates of x. For these elements 
a uniform inf-sup condition is satisfied (see [Jj), that is, there exists a constant 
P > independent of the mesh grid size h such that 

. , (q h ,V-v h ) 

inf sup - — — — > p. (11) 

The approximate velocity belongs to the discrete divergence-free space 

v h ,r = x h , r n { X h e ffo(0/») d I (gh. v • x/0 - o Vg* g Q h>r -i} , 

which is not a subspace of V. 

Let (u,p) G (H 2 (Q.) d DV)x (iJ x (0)/R) be the solution of a Stokes problem 
with right-hand side g, we will denote by s/j = S/,(it) G the so-called Stokes 
projection (see (20]) defined as the velocity component of the solution of the 
following problem: find {sh,<lh) £ {Xh, r , Qh,r-i) such that 

KVs/,, V<p h ) + (Vfffc, 0/0 = (g, 0/0 V</» h G X h , r , (12) 

(V-» fc ,Vfc)=° VV'ft G Qh,r-i- (13) 

The following bound holds for 2 < Z < r: 

\\u - fl fc ||o + fe||u - s h ||i < C/i'dluH, + ||p|Ui-i/b). (14) 

The proof of (fT4]l for = fi/j can be found in [3D]. The bound for the pressure 
is [IB] 

lb - ftllLVR ^ C/sfr' _1 (Nli + Ibiu.-VR), (is) 

where the constant Cp depends on the constant ft in the inf-sup condition pip . 



5 



We consider the semi-discrete finite-element approximation (uh,Ph) to (u,p), 
solution of (P)-©. That is, given uu (0) = HhUo, we compute ti/f(i) G Xn, r 
and pnit) G Qn,r-i, t G (0,T], satisfying 

(liff, </>if) + f(Vliff, V<^h) + b{u H ,U Hl 4> H ) + (Vp_y, = (/, (j) H ) y<f>H G -X"ff,r, 

(16) 

(V-u ff ,^) = V^h G Qh,t-1, 

(17) 

where 6(it, u, uj) = ((u ■ V)u + ^(V • u)v, w) for any u,v,w G Hq {Q) d - 

For 2 < r < 5, provided that (fT3]> -tfT5 ]) hold for I < r, and (|9]l-(|T0)) hold for 
k = r, then we have 

||«W-«H(*)||o + ^IKt)-« ff (*)l|i<C' i5 ^ 7 5, o<t<r, (is) 

(see, e.g., [H H HO]), and also, 

lb(t)-p«(*)IUvK^ C '^iV2' 0<t<T, (19) 
where r' = r if r < 4 and r 1 = r + 1 if r = 5. 

3 The new postprocessed method 

The postprocessing technique we propose is a two-level or two-grid method. In 
the first level, we choose a coarse mesh of size H and compute the mixed finite- 
element approximation (uh,Ph) to (u,p) defined by (I16ll - (fl~7|) . In the second 
level, the discrete velocity and pressure (uH(t),pii(t)) are postprocessed by 
solving the following linear Oseen problem: find (uh(t),ph(t)) £ {Xh,r, Qh,r-i), 
h < H, satisfying for all 4>h G Xh. r and iph G Qh,r-i 

v{Vu h {t),(f> h ) + ((u H (t) • V)S h (t), <t> h ) + (Vpfc(<), <t>h) = (/(*) ~ ( 2 °) 

(V-fihW.V'h) = 0- (21) 

Equations (|2"Uj) - (|2"Tj) can also be solved over a higher order mixed finite-element 
space over the same grid. For simplicity in the exposition we will only consider 
the case in which we refine the mesh at the postprocessing step. 

Let us observe that projecting equation (I2U1) over the discretely-free space 
Vh,ri and avoiding for simplicity the dependence on t in the notation, we get 
that Uh G Vh, r satisfies 

l/(VUh,Vh) + {{u H ■ V)u h ,v h ) = (/ - uhAh), VUfc G Vh,r- (22) 

We now prove that equation (|2"2"|) is well-posed, i.e., for H small enough there 
exists a unique function Uh G Vh,r solving (|2"2"|) . Let us denote by B H the bilinear 
form defined by 

B H (u h , v h ) = v{Vu h , Vv h ) + ((u H ■ V)u h , v h ), u h , v h G Vh, r - (23) 
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We proceed to show that B H is coercive which implies that there exists a unique 
function Uh € Vh,r satisfying (|22j) . Let us also observe that once a unique Uh 
is found, using the inf-sup condition flTlT) one easily obtains the existence and 
uniqueness of the pair (v,h,ph) satisfying (|2T)|) - (j2"l"j) . 

Lemma 1 Let Bh be the bilinear form defined in i23\) . Then, there exists a 
constant C such that for t > the following bound holds: 

\B H (v h ,v h )\ > U-C t(r _ 2)/2 )\\v h \\l Vv h GV h , r , (24) 

where 7 = 1/2 if the dimension d is d = 2, and 7 = 1/4 if d = 3. 

Proof To prove the coercivity of B H we follow [25j p. 2042]. Let us first 
observe that for any Vh £ Vh jr 

B H (v h ,v h ) = u\\Wvh\\o - ~(V • u H ,v h ■ Vh). 

Let qn be the I? orthogonal projection of Vf t ■ Vh over Qh,t-1) so that applying 
standard finite-element theory [5] and interpolation theory on Hilbert spaces (see 
e. g. ED § n -2] we have \\v h ■ v h - q H \\i^(p)/9. < CH~< \\v h ■ v h \\ 7 , for 7 G (0, 1]. 
Taking into account that the velocity u satisfies V • u — then 

B H (v h , v h ) = u\\Vv h \\l - -(V • (u H - u), v h ■ v h - q H ). (25) 

And then 

|(V • (u H - u),v h ■ v h - qu)\ < C\\u H - u\\i\\v h ■ v h - te||L2 (n)/R . 
Following [25, p. 2042] we get 

K'Vfc|| 7 <CK||?, (26) 

where 7 = 1/2 if d = 2, and 7 = 1/4 if d = 3. Using (gBJ together with ([18]) we 
get 

|(V • (u H - u),v h ■ v h - q H )\ < C^^inWvhWl 

Finally, going back to (f2"5"j) we reach (|2^|) . □ 

Let us observe that, for t > and H < (i( r ~ 2 ) '/ 2 v / 'C) 1 ' / ( r_1 +7) ; as a consequence 
of Lemma [D there exists a unique Uh £ V^.r satisfying (1221) . 

We introduce now a linearized problem that will be used in the proof of 
Theorem 1 where we state the rate of convergence of the new method. Let u 
be the velocity in the solution (u,p) of (d}-®. We will denote by (v,j) the 
solution of the following linearized problem 

- vAv + (u ■ V)v + S7j = d (27) 
div(u) = 
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in the domain ft subject to homogeneous Dirichlet boundary conditions. Let us 
observe that since the divergence of u is zero the bilinear form: 

B(v, w) — z/(Vu, Vw) + ((u ■ V)u, w), v, w G V. 

associated to this problem is continuous and coercive. Since the solution v G V 
of (|27|) satishes 

B(v,w) — (d,w), \/w G V 

by the Lax-Milgram theorem there exists a unique solution v. Due to ^ there 
exists also a unique pressure j. 

We will assume in the sequel that both problem ([27)) and its dual problem 
satisfy the regularity assumption 

IM|2 + ||j|| ff i ( fi)/ ffi <q|d||o. (28) 

The regularity assumption (|28|) can be proved by using the analogous regularity 
of the Stokes problem and a bootstrap argument, see [25j Remark 2.1]. 

In the following lemma we state the rate of convergence of the mixed finite- 
element approximation to the solution (v,j) of (|27|) defined as follows: find 
{vh,jh) G {X h ,r,Qh,r-l) such that 

v(Vv h ,V4> h ) + ((u-V)v h ,4> h ) + (Vj h ,<l> h ) = (d,^),V^6l V) (29) 

(V-v h ,i> h ) = 0, W> h G Q fc , r _i.(30) 



Lemma 2 Let (v,j) be the solution of and let (vh,jh) be its mixed finite- 
element approximation. Then, the following bounds hold for 2 < I < r 

||«-«fc||o + ^||«-«fc||i < Ch l (H|, + \\j\\m-*/n) , (31) 
\\j-3hhyn < Ch 1 - 1 (IMI; + ||i||H'-va) • ( 32 ) 

Proof Let us denote by Sh — Sh(v) the Stokes projection of v. More precisely, 
{sh,qh) G (Xh, r , Qh.r-i) will be the solution of (p~2|) - (fl3]) with right-hand-side 
g = d — (u ■ V)u. Let us denote by eh = Sh — Vh- Then, from (|29f and (jT2|) we 

get 

v(Ve h , VH);,) + ((« ■ V)e h ,w h ) = ((u- V)(s h - v),w h ), Vw^ G Vh.r- (33) 
Taking Wh — e^ in (|3"3"|) and using §5§ we get 

v\\eh\\\ < C\\u\\ L 2d/(d-i)\\s h - v\\i\\e h \\ L 2d < C||n||i/ 2 ||s/, - u||i||e/,||i, 
so that 

||efc||i<C||« h -t>||i. (34) 

Since \\v - v h \\i < \\v - s h \\i + \\eh\\i applying (jT4j) we conclude h\\v - v h \\i 
is bounded by the righ-hand side of (f3"Tj) . The bound (02) for the pressure is 
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readily obtained by means of the auxiliary value ku = Qh — jh- Subtracting (|29p 
from (|3Tj) and applying the inf-sup condition (fTTj) one easily gets 

/3||fcfclUVK ^ Hl e '*lli + C||w||i/2||sh - 

so that due to ([34]) and (jT4j) it follows that ||fc^||_L 2 /K is bounded by the right- 
hand side of ([211). Since \\j -jh\\v*/n < \\j ~ Qh\\ & /n + \\kh\\^/«., a PP!y m g C3 
we finally prove (J55]l. 

We are left with the task of proving the bound for the L 2 norm of the error 
in the velocity. We will argue by duality. Let us observe that 

II II \( e h,<P)\ for\ 

M [|o = sup m n — ■ (35) 
•pel 2 \m\o 

Let us fix ip G L 2 and let us denote by (w, k) the solution of the linearized dual 
problem 

—vAw — (u ■ V)ii! + Vfc = ip, 1 . „ 

div(iu) =0, J m ' (36) 

it = 0, on <9f2. 

As stated before we assume that this problem satisfies the regularity assump- 
tion ([55]). so that 

IM| 2 + ||fc||Hw K <qMlo- ( 37 ) 

We will denote by (wh,kh) £ {Xh t r,Qh,r-i) the mixed finite-element approxi- 
mations to (w,k). Reasoning exactly as before and applying ([37]) we obtain 

\\w-WhWx < Cfe(||Hla + ||A:|| H i/R) <Ch\\<p\\ , (38) 
||*-*fc||iVR < Ch (\\w\\ 2 + \\k\\ H1/R ) <Ch\\<p\\ Q . (39) 

Integrating by parts we reach 

(e h ,<p) = i/(Ve h ,Vii7) + ((u-V)e h ,«>)-((V-e h ),fc) 

= v(Ve h , V(to - to/,)) + ((u • V)e, i7 u) - t0 h ) - ((V • e h ),k- k h ) 
+v(Ve h , Vwh) + ((u ■ V)e h , w h ). 

And then, applying (j3"5|) and ([3"9"]) we reach 

|(e h) <p)\ < CHKI|i/>IMIo + C||«|| 1/2 [[e fc [|ift||^|| + CWe^ChMo 

+ \v(\7e h ,\7w h ) + ((« • V)e 7l , (40) 

Then, to conclude, it only remains to bound |^(Ve^, Vwh) + ((it • V)eft,,Wft)| 
which by ([33]) is equal to • V)(s^ — u), Let us decompose 

\((u ■ V)(s h - v),w h )\ < \((u ■ V)(s h - v),w h - w)\ + \((u ■ V)(s h - v),w)\. 

Then, integrating by parts in the last term 

\((u ■ V)(s h - v),w h )\ < C\\u\\ 1 / 2 \\s h - v\\i\\w h - H|i + \(( u • V)w,Sh - v)\, 
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and the bound for the first term on the right hand side above concludes by 
applying (fT4"f and (f38)l . Finally, since 



|((u • V)w, s h -v)\ < C\\u\\ L 2d/(d-i, \\Vw\\ L 2d\\s h - u|| . 
Applying Sobolev inequality ([5|) together with (|37[) and (fT4| we reach 

|((« • V)w, S h - V)\ < C||u|| 1/a |M| fc' + \\j\\ H l-i/K) , 

so that the proof is finished. □ 

We now state some results that will be use to get the rate of convergence of the 
new postprocessed method. The proof of the following lemma can be found in 
[TSJ Lemma 4] for the case r = 2 and in [TU1 Lemma 5.1] for r — 3,4. 

Lemma 3 Let (u,p) be the solution of {J)-© an d let u h be the mixed finite- 
element approximation to u. Then, there exists a positive constant C such that 

||«t(t)-M*)ll-i < ^^# r H°g(#)f\ te(0,T], r = 2,3,4 ) (41) 

HA^n («*(*)-«*(*)) ||o < ^^ff^llog^)!, te(0,T\, r = 3,4, (42) 

where r' = 2 when r = 2 and r' = 1 otherwise. 

The proof of the following lemma can be found in jTOl p. 226]. 

Lemma 4 Let (u,p) be the solution of ((J)-© and let ur be the mixed finite- 
element approximation to u. Then, there exists a positive constant C such that 

Kt)-tt H (t)||-i < : ^ m H T+1 \]ag{H)\,ts^,T\, r = 3,4. (43) 

We end this section with a theorem that states the rate of convergence of the 
new postprocessed method. 

Theorem 1 Let (u,p) be the solution of (H)-© and for r = 2,3,4 let ©-([10) 
hold with k = r + 2. Then, there exist a positive constant C such that the 
new postprocessed approximation (v,h(t) , Ph(t)) defined by i20 \) - ^l]) satisfies the 
following bounds for t £ (0, T] and H small enough: 

||u(t)-fi h (t)||i < Ch+-^H 2 \\og(H)\ 2 , r = 2, (44) 

\\u(t)-u h (t)\\j < j^h'-' + - ( ^ir+ 1 -i\]og(H)\, J=0,1, r = 3,445) 

||p(t)-^(t)|U VK < ^ 7 2^- 1 + ^ 7 2ff r |log( J H")r', r = 2,3,4, (46) 

where r' = 2 /or r = 2 and r' = 1 otherwise. 
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Proof Let us consider the linearized problem (|27p with right hand side d = 
f — tit. Then, the solution (v,j) of |27|) is the solution (u,p) of ([lj (J2j) . Let us 
denote by (i>h,jh) its mixed finite-element approximation, that is the solution 
of (I5TJ)) . This approximation satisfy the error bounds ([51]) and (J32]) for I = r. 
Let us decompose it-% = (tt - u^) + («/, - u^) andp-p/j = {p~jh) + {jh-ph)- 
To bound the first terms in these two decompositions we will apply (pTTj) and 
(|32|) . In the rest of the proof we deal with the other two terms. 

Let us denote by e/j = Vh — Uh- Subtracting (j2"0)) from (f2T)| it is easy to see 
that eh satisfies 

v(S7eh, V0h) + {{uh ■ V)eh,4>h) = {uh - u t , 4>h) + {{{uh - u) ■ V)v h ,4>h), 

for all 4>h G Vh,r- Taking (\>h = eh in the above equation and applying (|24[) we 
get that for H < (i('-2)/2 z/ / c -)i/(r-i+ 7 ) there existg a constant q such that 

1 1 e /i 1 1 1 < C(\\u t - u H \\-i + \\u H - u\\i\\v h - u\\i + \\u H - u\\o\\u\\ 3 / 2 ), 
and applying (j4"Tj) from Lemma [3j (Tl8|) and ([3T|) we get 

H^lli < WW* HT llosWr' + ^(^(Mr- 1 + /T), (47) 

from which (j44|) and the case j = 1 in (|45j) are concluded. 

We now get the error bound for the pressure. Let us denote rh — jh — Vh- 
Subtracting (|20|) from (|29|) and using (fTT|) it is easy to obtain 

^IMUvr ^ HKIIi + cII^h lli^lhilli + \\uh - uh\\-i 

+ C\\u H - u||i||Uft - u\\i + C\\u H - u||o||tt||3/2, 

from which we get flU applying (gTJ), flU]) from Lemma (0 and (f^Tj) . 

To conclude we get the error bound for the velocity in the L 2 norm. We 
will argue as in the proof of Lemma [21 that is, recalling ([35]) . for (p G L 2 (f2) 
we consider the solution {us, k) of (|36|) . so that (|40|) holds, and we are left to 
estimate |^(Ve^, Vui/,) + ((u • V)e/j,u>h)|. It is easy to see that 

v{Ve h , Vw h ) + ((it • V)e h , w h ) = {u H ~ u u w h ) + {{{u H - u) ■ V)u h ,w h ). (48) 

Let us now bound the two terms on the right hand side of (l48l) . For the first 
one, using (|37| and (|38| we get 

{uH-u t ,w h ) = (uh ~ut,Wh- w) + (u H - u t ,w) 

< \\UH - Ut\\-l\\w h - w||l + m _1 n(ltH - W t )||o||^M|o 

< C||u H - ut||_i/i|^||o + CWA-^uh - ut)||o|M|o- 

Applying now gl]) and ((42]) we have that (u H -u t , w h ) is 0{H r+1 1 \og{H)\/ {t^-^' 2 ) 
for r = 3, 4. Finally, we will bound the second term on the right hand side of 
(|48|) . To this end we decompose 

{{{u H - u) -V)u h ,w h ) = {{{u H - u) -V)(% - u),w h ) + {{{uh - u) -V)u,w h ) 

< C\\u H - u\\i_/ 2 \\u h - u||i||w^||i + {{{u H - u) ■ V)u,w h ) 

< CH\\u h - u\\i\\(p\\o + {{{uh - u) ■ V)u,w h ), 
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where in the last inequality we have applied (fT5|) and we have bounded ||n>/i||i < 
C||io|| i < C 1 1 V 3 1 1 o - Then, to conclude, it only remains to bound — u) ■ 

V)u,Wh)- Adding and subtracting w we get 

{((u H - u) ■ V)u, w h ) = {{{u H - u) -V)u,Wh - w) + (((u H - u) -V)u,w) 

< C\\u H - u|| ||u||3/2||^/i - u>||i + C\\u H - u||_i||Vu • If ||i 

< C\\uh - u||o/l||y||o + C\\uh - u||-l||u||2||H|2 

< C\\UH - u||o/l||y||o + C\\u H - U||-l||^||0) 

where we have applied {57]) . To conclude we apply (fT5)l and Lemma |U □ 

Remark 1 We observe from Theorem Q] that the postprocessed method in- 
creases the rate of convergence of the Galerkin method in one unit in terms 
of H, the size of the coarse mesh. In the case of linear elements the improve- 
ment is only achieved in the H 1 norm of the velocity but it is not obtained 
in the L 2 norm. Analogous results had been obtained for the standard post- 
processing in the linear case, see [3], [IS]. Let us also observe that a correct 
selection of the coarse and fine mesh diameters gives for the new postprocessed 
method the same rate of convergence than the Galerkin method over the fine 
mesh, although, of course, with different constants in the error bounds. The 
advantage of the method we propose is the saving in computational effort. For 
the method we propose the time integration is performed using the standard 
Galerkin method over the coarse mesh and only at the final time we solve one 
linearized Oseen-type problem over the fine mesh. Let us observe that, for ex- 
ample, the selection H = h 1 / 2 allows to get for the new postprocessed method 
the rate of convergence of the fine mesh in the H 1 norm when using linear ele- 
ments. The selection H = h 2 ' 3 allows to get the rate of convergence of the fine 
mesh in the H 1 norm when using quadratic elements, the choice H = /i 3 / 4 gives 
the rate of convergence of the fine mesh in the L 2 norm also for quadratics and 
so on. 

The reason why we have not carried out the error analysis for higher than 
cubic finite elements is that, as in the papers [20] and [10J . due to the loss of 
regularity at t = no better than 0(H 5 \ log(£T)|) error bounds can be proved. 



4 Numerical experiments 

We consider the Navier-Stokes equations ([1} in the domain f2 = [0, 1] x [0, 1] 
subject to homogeneous Dirichlet boundary conditions. For the numerical ex- 
periments of this section we approximate the equations using the mini-element 
[5] over a regular triangulation of fi induced by the set of nodes (i/N,j/N), 
< i,j < N, where N — 1/H is an integer. We study the spatially semi- 
discrete case. Hence, in the time integration (with the trapezoidal rule) suffi- 
ciently small time steps were taken so as to ensure that errors arising from the 
spatial discretization were dominant. In the first experiment we take the forcing 



12 



term f{t,x) such that the solution of (HJ-© with v = 0.05 is 

u (x,y,t) = irt sin 2 (7ra) sin(27ry), 
u 2 (x,y,t) = —Trtsm 2 (Try)sm(2irx), 
P(x,y,t) = 20tx 2 y. 

When using the mini-element it has been observed and reported in the literature 
(see for instance [32], [33], [5] [23], [28] and [29]) that the linear part of the 
approximation to the velocity, tA, is a better approximation to the solution u 
than Uh itself. The bubble part of the approximation is only introduced for 
stability reasons and does not improve the approximation to the velocity and 
pressure terms. For this reason in the numerical experiments of this section 
we only consider the errors in the linear approximation to the velocity. Also, 
following [3J, we postprocess only the linear approximation to the velocity, i.e., 
we solve problem (l2"0"l) - (|2"Tj) substituting ur and uh by u l H and u l H respectively. 
The finite element space at the postprocessed step is the same mini-element 
defined over a refined mesh of size h small enough to capture the asymptotic rate 
of convergence in the fine grid. The coarse and fine mesh sizes in the experiments 
arc H = 1/6, H = 1/8, H = 1/10 and H = 1/20 and h = 1/20, h = 1/26, 
h = 1/32 and h — 1/36 respectively. For the postprocessed approximation we 
also keep only the linear part. We apply the postprocessing step only once at 
time t = 0.5. In Figure [T] we have represented the size of the Galerkin and 




Figure 1: Galerkin (solid lines) and postprocessed (dashed lines) errors in L 2 
(asterisks) and H 1 (circles) for H = 1/6, 1/8, 1/10 and 1/12 and h = 1/20, 
1/26, 1/32 and 1/36 respectively. On the left, errors for the first component of 
the velocity. On the right, L 2 errors for the pressure. 

postprocessed errors with respect to the inverse of the coarse mesh size H. On 
the left part of the picture we present the results corresponding to the first 
component of the velocity. The results obtained for the second component of 
the velocity are analogous. On the right part of the picture we present the 
errors in the pressure. In both pictures, we have used solid line for the Galerkin 
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method and dashed line for the postprocessed method. The errors are measured 
in both the H 1 norm and the L 2 norm. In the picture, they are represented 
by circles (H 1 norm errors) and asterisks (L 2 norm errors). We can observe 
on the left of Figure [T] that, in agreement with the theory, the postprocessed 
method using the mini-element does not increase the rate of convergence in the 
L 2 norm of the velocity although the size of the errors are reduced. In the 
H 1 norm, however, also as predicted by the theory, the postprocessed method 
does increase the order of convergence by one unit (indeed, the errors of the 
postprocessed method in the H 1 norm are slightly smaller than those of the 
Galcrkin method in the L 2 norm). The same improvement is observed for the 
L 2 errors of the pressure on the right of Figure [TJ This means that we can 
obtain the level of error corresponding to the fine mesh at essentially the cost of 
the computation in the coarse mesh because the computation on the fine mesh 
is performed only once at time t = 0.5. Then, the dominant computational cost 
is caused by the time evolution in the coarse mesh saving time when compared 
with the time evolution in the fine mesh that is needed in a standard approach. 




Figure 2: First component of the velocity for the Galerkin method with v = 0.01 
and H = 1/10. 

In the next experiment we will show that the new postprocessed method 
produces better results than both the Galerkin method and the standard post- 
processed method ([4|. We consider now equations (Q} with initial condition 

u (x f y,t) = — 6 sin(7rx) 3 sin(7ry) 2 cos(7ry), 
u 2 (x,y,t) — 6 sin(7ra;) 2 sin(7r?/) 3 cos(7ra;), 

and forcing term / = 0. We take first v = 0.01. In Figured] we have represented 
the linear part of the first component of the velocity for the Galerkin method 
with H = 1/10 at time T — 0.5. In Figure[3]we show the standard postprocessed 
approximation with H = 1/10 and h — 1/30. We can observe that the standard 
postprocessing introduces some oscillations that were not present in the Galerkin 
approximation. These oscillations are not reduced with a smaller value of h. 
Finally, in Figure U we have represented the linear part of the first component 
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Figure 3: First component of the velocity for the postprocessed method with 
v = 0.01, H = 1/10 and h = 1/30. 



I 




Figure 4: First component of the velocity for the new postprocessed method 
with v = 0.01, H = 1/10 and h = 1/30. 

of the velocity for the new postprocessed approximation and the same values 
of coarse and fine mesh sizes, H = 1/10 and h = 1/30. We observe that this 
approximation does not oscillate at all and it improves the accuracy of both 
Galerkin and standard postprocessed approximations. 

In the last experiment we repeat the experiment with a smaller value of the 
diffusion parameter, v = 0.005, and the same values of H and h as that of 
the previous experiment. As it was already observed in the case of convection- 
diffusion equations [H] the behavior of the standard postprocessed method de- 
teriorates as the diffusion parameter decreases. 

We can observe that both the Galerkin and the new postprocessed approxi- 
mations, see Figures [5] and [7] respectively, do not present oscillations. As before, 
we can also observe the smoothing effect achieved by postprocessing with the 
new method proposed in this paper. On the other hand, the standard postpro- 
cessed method produces a completely wrong approximation, see FigureEl Let us 
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Figure 5: First component of the velocity for the Galerkin method with v = 
0.005 and H = 1/10. 




Figure 6: First component of the velocity for the postprocessed method with 
v = 0.005, H = 1/10 and h = 1/30. 

remark that, as it has been noted before in the literature, see [6], [30], the bubble 
functions used in the mini-element to generate a stable mixed finite-clement sat- 
isfying the inf-sup condition (|11[) have also a slightly stabilizing (over-diffusive) 
effect for moderate values of the Reynolds number. This fact explains the non- 
oscillating behavior of the linear part of the approximation to the velocity in 
the Galerkin method of Figure [5] We can also observe, see Figure [3 that the 
over- diffusive effect appearing in the Galerkin approximation of Figure O is at- 
tenuated by postprocessing with the new method. 
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Figure 7: First component of the velocity for the new postprocessed method 
with v = 0.005, H = 1/10 and h = 1/30. 
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